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Abstract 

We analyze the stability of the Einstein static closed and open universe in two types of exponential 
f(T) gravity theories. We show that the stable solutions exist in these two models. In particular, 
we find that large regions of parameter space in equation of state w = p/p for the stable universe 
are allowed in the f(T) theories. 
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I. INTRODUCTION 



The Einstein static universe has recently been revived since our universe might evolve 
from it to inflation. This is called the emergent universe with an inflationary singularity to 
avoid a big bang singularity 

One way to look at the theory beyond general relativity (GR) is the teleparallel equivalence 



n 



of general relativity (TEGR) |2|-|6| introduced first by Albert Einstein pj . Different from GR 
with the Levi-Civita connection, teleparallel gravity (TG) uses the Weitzenbock connection, 
which has no curvature (R) but has torsion (T). To explain the late-time acceleration of 
the universe, a modified gravity has been proposed by extending T in the TG action to an 
arbitrary function f(T). In f(T) gravity, torsion is responsible for the acceleration of the 
universe {si, 9]. This modified gravity would not only avoid the big bang singularity 
but also provide an alternative to inflation ll|. However, /( T) gravity has some intrinsic 
problems, such as the violation of local Lorentz invariance 12[. In some of modified gravity 
models, the solutions of the stable Einstein static universe do exist under linear homogeneous 

. In Ref . [15j , the stability of the Einstein static 
universe in a power law f(T) plus the cosmo logical constant A is studied. However, A is at 
the pre-inflation scale, which may lead to the hierarchy problem in the current universe. To 
avoid this problem, we examine the stability of the Einstein static universe in two kinds of 
exponential f(T) gravity theories {{J, which could explain the static Einstein universe in the 
pre-inflation era and be reduced back to TEGR in the present universe. 

The paper is organized as follows. In Sec. [Til we briefly introduce the formulation of f(T) 
gravity with the Weitzenbock connection. We also show two explicit forms of the exponential 
theories. In Sec. IHIt we give the conditions for the stable Einstein static solutions. In Sec. IIVI 
we study the stable solutions in the exponential models. We illustrate the oscillating behavior 
of the Einstein universe in Sec. [V] We present our conclusions in Sec. IVII 



II. FORMULATION IN f(T) GRAVITY 

In teleparallelism, the dynamical object is the vierbein field ej(x M ), which is an orthogonal 
basis for the tangent space at the point x M of the manifold with the relation: e, ■ e 3 - = r]ij, 
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where 77^ =diag(l,- 1,-1,-1). The metric tensor is given by 

9n» = Vij^( x ) e l( x ) (2- 1 ) 

with 

ei(x)cf(x) = 5) . (2.2) 
In this formulation, the Weitzenbock connection is used and the torsion tensor is defined by 

Z£ = " r £, = ^{9 A - 9 v el) , (2.3) 

where 



The action in TG is expressed as 



rj, = e&ej, . (2.4) 
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where e = det(e t ) = y/—g and T is the torsion scalar, defined by 



T = S x ^T x flu , (2.6) 

with 

V" = \{K"\ + 5\T e \ - S\T e %) , (2.7) 

and 

K<+\ = ^j-(T^\-T u \-T x ^). (2.8) 
The modified teleparallel action for /(T) gravity is given by 9] 

,= ikS ef(T)i,x ' (2 - 9) 

where f(T) is an arbitrary function of T. In this paper, we will concentrate on the following 
two types of the exponential f(T) models: 



f(T) = T + aT (1 - e^ To/T ) , (2.10) 

and 

/(T) = T + «T (l-e^o), (2.11) 

where a, (5 and T are constants. The motivation of these two models is that f(T) can be 
reduced back to TEGR for small T, corresponding to the current universe. Explicitly, when 
T is very small compared to To and < 0, the two models in Eqs. (12.101) and (12. lip give 
}'{T) w (1 + a)T and f(T) ~ T, respectively. 



III. EINSTEIN STATIC UNIVERSE 



We now consider the FRW metric 



d s 2 = dt 2 - k 2 a 2 (t)[d(kip) 2 + sm\kip)(d9 2 + sin 2 9d(j) 2 )] 



(3.1) 



where its vierbein fields are given by 16] 



p- —i p° — p° — p u — p L — p 2 — p 6 — n 

e o ~~ c ip ~ c e ~ c </> ~~ e o ~~ e o ~~ e o ~~ u ' 

e]p =—a(t)k 2 cos9, el = a(t)k sm(kip) sin 9 cos (kip), el = — a(t)/csin 2 (/c^) sin 2 

e 2 =a{t)k 2 sin # cos 0, e 2 , = —a{t)k sin(fc-0) [sin(/r0) sin — cos{kip) cos 6 1 cos 0] , 

e 2 =— a(£)fc sin(/r0) sin # [cos^VO s i n + sin(fc^) cos 9 cos 0] , 

e| =—a[t)k 2 sin sin 0, = — a(t)A; sin(£;0>) [sia(A;'0) cos + cos(Aj^) cos 9 sin 0] , 

e\ =—a(t)ksm(kip) sin # [cos(kip) cos — sin(A;^) cos 9 sin 0] , (3-2) 



with k = 1 and z or /c 2 = ±1, representing the closed and open universe, respectively. The 
torsion scalar can be written as 



T = Q(k 2 a~ 2 - H 



As a result, the modified Friedmann equation is given by 16] 



(3.3) 



12 Hf'(T) + f(T) = IQnGp = K p, 



(3.4) 



(k 2 a~ 2 + H)(48H 2 f"(T) + Af'(T)) - f'(T)(8H + 12H 2 ) - f(T) = k P . 



(3.5) 



To get an Einstein static universe, the conditions of a = H = 0, a = 0, and To = T(ao) 
fc 2 6/aQ are imposed. By using Eqs. (13. 4p and (13. 5p . we obtain 



fo = f( T o) = K Po 

k —5" - JO = Kpo 



(3.6) 
(3.7) 



with /q = df /dT\ T=To , p = p(a ) and p = p(°o)- By combining Eqs. (13. 6 p and (13. 7p . we 
find 

'Tf r 



= |(l+w), 

J J T=T Z 

where w = p/p is equation of state of the background matter. 



(3.8) 



We now perform the linear homogeneous scalar perturbations in the static Einstein uni- 
verse. The perturbations in a and p depend only on time, i.e., 



a(t) = a (l + Sa{t)), p(t) = p Q (l + 6p(t)) 



with 



15| 



5a(t) = de 7 ' + C 2 e-^ 
By using the procedure in Ref. jl5|, we get 



2 K P0 ft , \ 

7 = ttttt(1 + w) 



4/<5 s 



;i + 3w)^-3(l + w; 



/o/o 

/J 



J 



(3.9) 



(3.10) 



(3.11) 



From Eq. (I3.10p . we obtain an oscillating universe if 7 2 < 0, corresponding to the stable 
Einstein static universe. 

In the following section, we discuss the stable Einstein static solutions of the exponential 
f(T) models in Eqs. (I2TTUD and (123T]) . 

IV. CONDITIONS FOR STABLE EINSTEIN STATIC SOLUTIONS 
A. The first type of the exponential models 

In the closed universe, we have To = Q% 2 for the Einstein static solution. Substituting 
Eq. fl2.10p into Eqs. (I3.6P and (13. 7ft . we obtain the constraints for the Einstein static universe 



To(l + a — ae^) = K,p 
^Toafie 13 = Kp f w + - 



(4.1) 
(4.2) 



leading to 



Kp 



1 - e 13 

1 - -r-^- (1 + 3w) > 



2(3e 



13 



(4.3) 



where T /(Kp ) is positive since po is positive definite and T > in the closed universe. 
From Eqs. fl2~T0l . fl3~TT|) . (Q]l and we have 



^(l + w)(l + 3w) 



T 



(4.4) 



Note that f(T) — > T + a (3 To when /3 — >• 0. If the combination a (3 To is finite, we may define 
it as the cosmological constant A. As expected, in the limit of f(T) — > T + A, one has 
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np 



l + w)(l + 3w), A 



3w) 



(4.5) 



which is the result in GR. For the stable universe in GR, which requires 7 < 0, one finds 

-l<w<--. (4.6) 

In Fig. [I] (left panel), we show the stable regions (7 2 < 0) based on Eqs. (14.31) and (I4.4p . For 
f3 = 0, as seen from the figure, we obtain the same allowed region of— 1 <w< —1/3 as that 
in Eq. (Oil . 

In the case of the open universe, from Eqs. (13.61) and (13. 7\) we get the constraints for the 
Einstein open static universe as: 

r (l + a - ae p ) = k Po , (4.7) 
2 T af3eP = k Pg ( w + \ ) , (4.8) 



3 ' V 3, 

respectively, where T = — 6oq 2 . With these two constraints, we find 

which is negative because T = — 6oq 2 < in the static open universe. From f{T) in 
Eq. (13. lip , we derive 



7 2 

' ST, 



^(l + w)(l + 3w)(l + ^ + w^) . (4.10) 



In Fig. [T] (right panel), we display the stable solutions (7 2 < 0) for the Einstein static open 
universe. Clearly, there is no stable region at (3 = 0. Thus, the universe is not stable in the 
spatially open universe. 

B. The second Type of the exponential models 

For the second type of the exponential models in Eq. (12. lip , we obtain conditions for the 
closed Einstein static universe as 

/o = T (l + a - ae?) = k Po , (4.11) 

|r (l - 2«/3e /3 ) - k Po = k Po , (4.12) 



leading to 



r _ |(l + w )(l-e0 + 2^ 



k Po 1 - + 2(3e? 
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FIG. 1. Stable solutions (shaded regions) of the first exponential gravity model in the spatially 
closed (left panel) and open (right panel) Einstein static universe. 

which will be positive because T = 6a^ 2 > in the closed universe. Substituting Eq. (12. lip 
into Eq. (13. lip , we get 



3k p 

1 = 5W? (1 + W) 



; 1 + w) (1 + 3w)- 2(1 + 2 ,)(w + 1 - < 1+ ^-^7 

114) 

From Eqs. (I4.13P and (I4.14p . we can plot the stable solutions (7 2 < 0) of the Einstein static 
closed universe in Fig. [2] (left panel) . 

In the open universe, the forms of equations are the same as those in the closed universe. 
However, since T = — 6ag 2 < in the static open universe, the constraint in Eq. (I4.15P is 
negative, i.e. 



< 0, (4.15) 



To_ _ f(l + w)(l-eP) + 2/^ 
~Kp Q ~ 1 - ef 5 + 2(3e? 

while 7 2 is also negative for the stable static universe: 

„ ^Vo , N \, ^ , o..a nn . ( (1 + w)(l - eP) + \^ 



7 = 77^(1 + W 

1 8T /o 21 



1 + w)(l + 3w) - 2(1 + 2/3) w + 1 



1 - c 3 + 2/3^ 

The stable solutions for the open universe are illustrated in Fig. [2] (right panel). 



.16) 



* o 





FIG. 2. Legend is the same as Fig. Q]but for the second model. 



V. OSCILLATING EINSTEIN UNIVERSE 



From Eq. (I3.10p . it is obvious that the stable Einstein static universe is oscillating. The 
behavior has also been explicitly demonstrated in other gravity theories, such as the vacuum 
energy model 1?| and DGP braneworld scenario [iuj . In this section, we would like to give 
similar discussions to illustrate this oscillating behavior in the first type of the exponential 
models in Eq. (I2.10p . The results can be easily extended to the second one. 

Since an obvious stable solution for the first exponential closed universe comes from (3 = 
—2 and w = 0, we can get 

from Eq. (14.41) . Consequently, the scale factor is given by 



a(t) = a 



1 + C sin 



8T /o 2 v/3 



+ rj 



(5.2) 



By defining t' 



&t, we get 



a(t') 
d(0 



a C 



1 + C sin I —= + 77 

V3 



( l ' 
_ cos — = 



V 



(5.3) 




Using initial conditions given by a(t = 0) = 1.1 x ao and a(t = 0) = 0, we obtain r] = n/2 
and C = 0.1. We plot the stable small oscillating closed universe for the first exponential 
gravity model in Fig. EJ 

For the open universe, we pick the stable solution of (5 = —4 and w = —1/2. With the 
same initial conditions as the closed universe, we find 



q(t') 
a 

a(t') 
a 



1 H sin 

10 



1 

10 




cos 




(5.4) 



where t' 
in GR. 



m\fo 



Tit. In Fig. HI we show the oscillation in the open case, which is not allowed 
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VI. CONCLUSIONS 



We have discussed the linear homogeneous scalar perturbations near the Einstein static 
universe in f(T) gravity. We have explicitly studied both closed and open universe in the two 
exponential /(T) gravity models. These two f(T) models are proposed in order to explain 
the pre-inflation universe and be reduced back to TEGR in the present universe. We have 
demonstrated that in these models, the Einstein static universe can be stable in both open 
and closed cases with large allowed regions of w. Note that in GR, only the closed Einstein 
universe contains stable solutions with — 1 < w < — ~ for the background matter. Explicitly, 
for the first exponential gravity model of f(T) = T + aT(l — e l3T °^ T ) 1 we have shown that 
w > — 1 (w < — |) in the closed (open) universe can have the stable solution. In the limit 
j3 — > 0, the model gives the GR result as f(T) — > T+A with A = a(3T being the cosmological 
constant. For the second exponential gravity model of f(T) = T + oTq(1 — e^ T2 / T o ), we have 
found that w > — 1 (arbitrary w) is the stable solution for the open (closed) Einstein universe. 
We have also illustrated the oscillating behaviors of the stable Einstein static universe in the 
first model. We have shown that both a(t) and a(t) are oscillating without divergences as 
expected. 
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